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ABSTRACT 

Quantum mechanical tunneling inversion transition in ammonia NH3 is actively used as a 
sensitive tool to study possible variations of the electron-to-proton mass ratio, fi — m c /m p . The 
molecule H 3 + has the inversion barrier significantly lower than that of NH 3 . Consequently, its 
tunneling transition occurs in the far-infrared (FIR) region and mixes with rotational transitions. 
Several such FIR and submillimiter transitions are observed from the interstellar medium in the 
Milky Way and in nearby galaxies. We show that the rest-frame frequencies of these transitions 
are very sensitive to the variation of /z, and that their sensitivity coefficients have different signs. 
Thus, H30 + can be used as an independent target to test hypothetical changes in [i measured at 
different ambient conditions of high (terrestrial) and low (interstellar medium) matter densities. 
The environmental dependence of \x and coupling constants is suggested in a class of chameleon- 
type scalar field models — candidates to dark energy carrier. 



Subject headings: molecular data 
elementary particles 

1. Introduction 



techniques: radial velocities — ISM: molecules — dark energy 



The spatial and temporal variability of dimen- 
sionless physical constants has become a topic of 
considerable interest in laboratory and astrophys- 
ical studies as a test of the Einstein equivalence 
principle of local position invariance (LPI), which 
states that outcomes of nongravitational experi- 
ments should be independent of their position in 
space-time (e.g., Dent 2008). The violation of 
LPI is anticipated in some extensions of Stan- 
dard Model and, in particular, in those dealing 
with dark energy (e.g., Hui et al. 2009; Damour & 



Donoghue 2010). 

A concept of dark energy with negative pres- 
sure (p = — p) appeared in physics long before 
the discovery of the accelerating universe through 
observations of nearby a nd distant (at redshif z ~ 
1) supernoyae ty pe la ( Perlmutter et al. 19981: 



Riess et al. I H998). Examples of dark energy in a 



form of a scalar field with a self- interaction poten- 



tial ca n be found in revi e ws bv lPeebles fc Ratra 



2003 ). ICopeland et al. I ([20061) . and by lUzan 



20101 ). Since that time many sophisticated mod- 



1 On leave from A. F. Ioffe Physical- Technical Institute, 
Saint Petersburg 194021, Russia 



els have been suggested to explain the nature of 
dark energy and among them the scalar fields 
which are ultra-light in cosmic vacuum but pos- 
sess a large mass locally when they are coupled to 



1 



ordinar y matter by the so-call e d chameleon mech- 



anism (IKhourv & Weltman 2004; Br ax et al 



2004tlAvelino Il2008t iBrax et al. Il2010allbh. A sub 



class o f these models considered bv lOlive fc Pospelov 



(|2008l ) predicts that fundamental physical quanti- 
ties such as elementary particle masses and low- 
energy coupling constants may also depend on the 
local matter density. Since the mass of the elec- 
tron m e is proportional to the Higgs vacuum ex- 
pectation value (VEV ~ 200 GeV), and the mass 
of the proton m p is proportional to the quan- 
tum chromodynamics (QCD) scale Aqcd ~ 220 
MeV, we may probe the ratio of the electroweak 
scale to the strong scale through the measure- 
ments of the dimensionless mass ratio /i = m e /m p 
in high density laboratory (terrestrial) environ- 
ment, jttiab, and in low density interstellar clouds, 
/"space (/Oiab//O sp aco > 10 10 ). In this way we 
are testing whether the scalar field models have 
chameleon-type interaction with ordinary mat- 
ter. Several possibilities to detect chameleons 
from astronomical observations were discussed 



where c is the speed of light and AV = V Iot — V lnv . 

Surprisingly, recent observations of a sample 
of nearby (distance R ~ 140 pc) cold molecu- 
lar cores (T kin ~ 10K, n = 10 4 - 10 5 cm" 3 , 
B < 10 fiG) in lines of NH 3 ( J,K) = (1,1) at 
24 GHz, HC 3 N J = 2 - 1 at 18 GHz, and N 2 H+ 
J = 1 — at 93 GHz reveal a statistically sig- 
nificant positive velocity offset between the low- 
J rotational and inversion transitions: AV = 



V Iot 



V im = 27 ± 4 stat ± 3 sys m s 1 , which gives 



(26±4 stat ±3 sys )xl0 9 (jLevshakov et al 



2010al i 1 1. A few molecular cores from this sample 



were mapped in the NH 3 (1,1) and HC 3 N (2-1) 
lines and it was found that in two of them (L1498 
and L1512) these lines trace the same material and 
show the offset AV = 26.9 ± 1.2 sta , t , ± 3.0 SVR m s" 1 



throug hout the entire clouds (|Levshakov et al 
12010b ). It was also demonstrated that for these 



in 


Burraere et al. 


(2009^ 


Brax & Zioutas 1 ( 


201(A. 



Davis et al. (2009) 



and lAvgoustidis et al. 



First experiments constraining these 
models were recently carried out in Fermilab 



mod els were recently carriecl out m fermiiab 
( Upadhve et al. Il2010f) and in Lawrence L ivermore 



2010h . 



National Laboratory (jRvbka et al. 

At the moment, the most accurate relative 
changes in the mass ratio A/i//i = (/J. spa cc — 

Mlab ) / Mi ab can be obtained with the ammonia 

meth od ( van Veldhoven et al. l2004lFlambaum fc Kozlov 
2007h . NH 3 is a molecule whose inversion fre- 
quencies are very sensitive to any changes in fi 
because of the quantum mechanical tunneling of 
the N atom through the plane of the H atoms. 
The sensitivity coefficient to //-variation of the 
NH 3 (J,K) — (1,1) inversion transition at 24 
GHz is Qi n v = 4.46. This means that the inver- 
sion frequency scales as Alu/lj = 4.46(A/i//i). In 
other words, sensitivity to ^-variation is 4.46 times 
higher than that of molecular rotational transi- 
tions, where Qrot = 1< Thus, by comparing the 
observed radial velocity of the inversion transition 
of NH 3 , Vl n v, with a suitable rotational transition, 
Vrot , of another molecule arising co-spatially with 
ammonia, a limit on the spatial variation of /i can 
be determined: 



clouds the frequency shifts caused by external elec- 
tric and magnetic fields and by the cosmic black 
body radiation-induced Stark effect are less than 1 
m s _1 . Optical depth effects in these clouds were 
studied from the analysis of unsaturated (r < 1) 
and slightly saturated (r « 1 — 2) hyperfine com- 
ponents of the corresponding molecular transitions 
and it was found that both groups of lines have 
similar radial velocities within the la uncertainty 
intervals. 

The nonzero Afi implies that at deep interstel- 
lar vacuum the electron-to-proton mass ratio in- 
3reases by ~ 3 x 10" 8 as compared with its terres- 
trial value and, hence, LPI is broken. In view of 
the potentially important application of this dis- 
crepancy to the fundamental physics, one has to 
be sure that the nonzero A/i is not caused by some 
overlooked systematic errors. An obvious way to 
tackle this problem is to use other molecular tran- 
sitions which have different sensitivity coefficients 
Qp. It has already been suggested to measure 
A-double t lines of ligh t diatomic molecules OH 
and CH ( Kozlov 1 12009^ . and microwave inversion- 
rotational transitions i n partly deuterated ammo- 



nia NH 2 D and ND 2 H (|Koz1ov et al. Il2010h 



In the present paper, we propose to use tun- 
neling and rotation transitions in the hydro- 
nium ion H 3 + . Like ammonia, it also has 
a double minimum vibrational potential. The 



A^ 



Vrot — Mnv 
c(Qhw Qrot) 



Ay 1 Presente d are the corrected values of AV and A/i//^ dis- 

0.3 , (1) cussed in lLevshakov et al. I ll2010bh . 
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Fig. 1. — The level scheme for H3O 
picted frequencies are in GHz. 



The de- 



inversion transitions occurs when the oxygen 
atom tunnels through the plane of the hydro- 
gen atoms. This leads to an inversion split- 
ting of the rotational levels. The splitting of 
H,cjO+ is very l arge, 55.3462 ± 0.0055 cm" 1 
(iLiu fc Oka l f"l985) as compared to 1.3 cm 1 split- 
ting in NH 3 . Consequently, the ground-state 
inversion-rotational spectrum of H3O" 1 " is observed 



in th e submillimeter-wave regio n (Plummcr et al. 



19851 iBogev et all 119851 IVerhoeve et al. Ill988l ). 



whereas pure i nversion 



infrar ed region (jVerhoeve et al 
20091) . 



transiti ons i n the far- 
19891 



Yu et al. 



H30 + has both ortho- and para-modifications 
(see Fig. [IJ. In the submillimeter, - the range 
accessible from high altitude ground-based tele- 
scopes, - there are three low-lying transitions at 
307, 364, and 388 GHz which belong to para- 
H30 + , and one ortho-H30 + transition at 396 
GHz. The 388 GHz line is, however, blocked 
by water vapor in the atmosphere. The other 
lines were observed i n the interstellar molecular 
1981 



clouds (Hollis et al. 



Woottcn et al 
Phillips et al 



1991 
1992h 



Wootten et al. 



1986; 



van der Tak et al. I 12006 ; 
The 364 GHz line was 
also observed in two galaxies: M 82 and Arp 220 



(|van der Tak et al. 1120081) . In far-IR, H 3 0+ lines 
were detected fro m aboard the space observa tories 
at u = 4.31 THz (iTimmermann et al. Ill996l). u 



1.66, 2.97, and 2.98 THz dGoicoechea fc Cernicharo 



200ll;lLerate et alj2006UPolehampton et al. II2007I) . 



984 GHz dGerin et al. I 120101) and 



1.03,1.07, and 1.63 THz (|Benz et al. lEoioh . 



The observed transitions of H30 + arise in the 
warm (T kin ~ 100 K) and dense (n « 10 5 - 
10 6 cm -3 ) star- forming regions surrounding pro- 
tostars where hydronium appears to be one of the 
most abundant species with the abundance as high 



as XfH3Q+) 5 x 10" 9 (jWootten et all Il99lt 



Benz et al. Il2010h 



2. Sensitivity coefficient of inversion tran- 
sition 

Sensitivity of the H30 + inversion transition 
to /i-variation can be estimated from the analyt- 
ical Wentzel-Kra mers-Brillouin (WKB) app roxi- 
mation. Following Landau fc Lifshitzl (|l977f ). we 
write for the inversion frequency (used units are 
H = |e| = m e = 1): 



2E, 



(2) 



where S is the action over classically forbidden 
region and Eq is the ground state vibrational en- 
ergy. Expression ^ gives the following sensitivity 
to /i-variation: 



S+l 



SE 



2(C/max — ^0) 



(3) 



where U max is the barrier hight and we are not 
using an additional approximation Eq — uj v /2. 

Acc ording to lRaiamaki et al. ( 2004 ) and lDong fc Nesbitt 
(|20061 ). we can take U max = 651 cm -1 and 
Eq « 400 cm -1 . The inversion frequency for 
H 3 0+ is 55.3 cm- 1 . Thus, Eqs. (2© give: 



1.5, 



2.5. 



(4) 



Dong fc Nesbitt I (|2006l ) report the inversion 



frequencies for H 3 0+, H 2 DO+, HD 2 0+, and 
D 3 0+ to be 55.3 cm" 1 , 40.5 cm" 1 , 27.0 cur 1 , 
and 15.4 cm -1 , respectively. Neglecting the weak 
dependence of the reduced mass, m r , on the in- 
version coordinate, we ta ke m r to be 0.7, 0.8, 1.0 , 
and 1.25, respectively (|Dong fc Nesbitt I 120061) . 
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Figure [2] shows the inversion frequency as a func- 
tion of m r . From this plot we can estimate the 
sensitivity coefficient for H30 + to be: 



2.46. 



(5) 



which is in a perfect agreement with Eq. We 
can conclude that the inversion transition in H30 + 
is almost two times less sensitive to /i-variation, 
than similar transition in N H 3 , where Qi nv = 4.5 
(iFlambaum fc Kozlov 112007ft . 



3. Sensitivity coefficients of mixed transi- 
tions 

The spectrum of the rotatio nal and inversio n 
transitions of H 3 + is studied in Yu et all ([2009). 
For the lowest vibrational state we can write the 
simplified inversion-rotational Hamiltonian as: 



H = BJ{J + 1) + (C - B)K 2 - Dj[J(J + l)f 
-D JK J(J + l)K 2 - D K K 4 + ... 
+ ^ {W a + WjJ(J + 1) + W K K 2 + ...} 



(6) 



Here we neglected higher terms of expansion in J 
and K; s = ±1 for symmetric and antisymmetric 
inversion state; total p arity p = ( — 1) K s. Numeri- 
cal values are given in lYu et al. I (|2009h (MHz): 



Note that we write Hamiltonian ([6]) in such a 
way that terms which determine inversion split- 
ting are collected in the last line. Therefore, we 
h ave the following relation with parameters used 
in 



Yu et al. (200 



B = [B(0+) + B(0-)]/2. 
Wj = S(0+) - B(0~) , 



(8) 



and similarly for C — B and Wk ■ Parameters Dj, 
Djk, and Dk are averaged over inversion states 
s = ±1. 

To estimate sensitivities of the mixed transi- 
tions it is sufficient to account for /i-dependence 
of the dominant parameters B, C, and Wq. It is 
clear that B, C ~ /i and Wq scales as yfi™ . It 
follows, that for rotational part of the energy we 




0.6 0.7 0.8 0.9 1.0 U 1.2 1.3 

Reduced mass (arb. units) 

Fig. 2. — Inversion frequency as a function of the 
reduced mass for hydronium ion isotopologues. 



have Qrot = 1 and for inversion part Q mv is given 
by Eq. (HJ) or (J3J) - Thi s leads to the expre ssions, 



used earlier for NH2D ( Kozlov et al. 2010|) 



B 


C-B 


Dj 


Dj K 






334406 


-148804 


35 


-70 


(7) 


and 


D K 


Wo 


Wj 


W K 




41 


-1659350 


5988 


-8458 







W ro t i Wj, 



W ro t 



± 



(9) 



(10) 



We use Hamiltonian © and expression (ITU1) to 
calculate the frequencies and sensitivities of the 
mixed transitions. The obtained results are pre- 
sented in Table [TJ Final results are very sensitive 
to the parameter Q lnv . A good agreement between 
two different estimates of Qi nv from Eq. (0| and 
from Fig. Q] shows that this parameter is known 
with 10% accuracy, or better. In the next ap- 
proximation, we need to weight independently all 
terms of the Hamiltonian © with different scal- 
ings. However, this does not lead to any significant 
changes in sensitivities of the low- J transitions 
from Table [U 

4. Discussion and conclusions 

We have shown above that the rest-frame fre- 
quencies of the inversion-rotational transitions of 
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Table 1: Frequencies and sensitivi ties to ji- variat i on of t he inversion-rota tion transitions in H 3 + . Experi 



mental frequencies are taken from lPickett et al. (1998); lYu et al. (2009) 







Transition 






Frequency (MHz) 




J 


K 


s 


J' 


K' 


s' 


Exper. 


error 


Eq. @ 




1 


1 


-1 


2 


1 


4-1 


307192.410 


0.05 


307072 


4-9.0 


3 


2 


+1 


2 


2 


-1 


364797.427 


0.10 


365046 


-5.7 


3 


1 


+1 


2 


1 


-1 


388458.641 


0.08 


389160 


-5.2 


3 





+1 


2 





-1 


396272.412 


0.06 


397198 


-5.1 








-1 


1 





+1 


984711.907 


0.30 


984690 


4-3.5 


4 


3 


+1 


3 


3 


-1 


1031293.738 


0.30 


1031664 


-1.4 


4 


2 


+1 


3 


2 


-1 


1069826.632 


0.30 


1071154 


-1.2 


3 


2 


-1 


3 


2 


4-1 


1621738.993 


2.00 


1621326 


4-2.5 


2 


1 


-1 


2 


1 


4-1 


1632090.98 




1631880 


4-2.5 


1 


1 


-1 


1 


1 


4-1 


1655833.910 


1.50 


1655832 


4-2.5 



H 3 + are very sensitive to the value of \x. For 
a given transition from Table [TJ Wj, with the 
sensitivity coefficient Qi, the expected frequency 
shift, Aujj/ujii due to a change in [i is given by 



( Levshakov et al. I ^OlOa) 



Au/j _ ojj - Uj _ A/i 

— — tyi 



(11) 



where Wi and Qi are the frequencies corresponding 
to the laboratory value of fj, and to an altered [i in 
a low-density environment, respectively. 

By analogy with Eq. ([T]), we can estimate the 
value of A/1//1 from two transitions with different 
sensitivity coefficients Qi and Qf 



A/i 



Vj-Vj 
c(Qi - Qj) ' 



(12) 



where Vj and Vi are the apparent radial velocities 
of the corresponding H 3 + transitions. 

Consider two lowest frequency transitions from 
Tableffl 17/ ->■ 2+ and 3^ ->■ 2~ of para-H 3 0+ 
at, respectively, 307 and 364 GHz. Here AQ = 
Q307 — Q364 — 14.7, which is 4 times larger then 
the AQ value from the ammonia method. This 
means that the offset AV^ ~ 27ms~ 1 , detected in 
the ammonia method, should correspond to the 
relative velocity shift between these transitions, 
AV = V/ 364 - V307, of about 100 m s _1 . 

Published results on interstellar H 3 + allow us 
to put an upper limit on A^if/j,. The observations 
of the 307, 364, and 396 GHz lines carried out at 
the 10.4-m telescope of the Caltech Submillimeter 



Observatory (CSO) by IPhillins et ail (|l992t) and 
the observations of the 307 and 364 GHz lines at 
the 12-m APEX telescope (Atacam a Path finder 
Experiment) by Ivan der Tak et al. (2006) have 
accuracy of about 1 km s _1 that provides a limit 
on A/x//i < 2 x 10 -7 , which is consistent with the 
signal ~ 3 x 10~ 8 revealed by the ammonia method 
((Levshakov et al. Il2010alibh . 

In order to check ammonia results we need to 
improve the accuracy of H 3 + observations by 
more than one order of magnitude. According 
to Table Q] the uncertainties of the laboratory 
frequencies of the transitions at 307, 364, and 
396 GHz are, respectively, 50, 80, and 45 m s _1 . 
Therefore, we also need a factor of few improve- 
ment of the laboratory accuracy to be able to de- 
tect reliably an expected signal AV ~ 100 m s _1 
and to check the non-zero ammonia results. 

An important advantage of the hydronium 
method is that it is based on only one molecule. In 
the ammonia method there is unavoidable Doppler 
noise caused by relative velocity shifts due to 
spatial segregation of NH 3 and other molecules. 
When using hydronium, the only source of the 
Doppler noise may arise from possible kinetic tem- 
perature fluctuations within the molecular cloud 
since the submillimeter H 3 + transitions have 
different upper level energies: E u = 80, 139, and 
169 K for the 307, 364, and 396 GHz transitions, 
respectively. It is also important that two H 3 + 
transitions have similar Q values: (5364 ~ £?396- 
This allows us to control the Doppler noise and 
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to measure accurately the relative position of the 
307 GHz line. 

The analysis of other possible sytematic effects 
for hydronium is mostly similar to what was done 
in detail for ammonia in Levshakov et al. (2010b). 
The systematic shifts caused by pressure effects 
are about a few m s^ 1 , or lower. As mentioned 
above in Sect. [TJ the frequency shifts caused by 
external electric and magnetic fields and by the 
cosmic black body radiation-induced Stark effect 
are less or about 1 m s _1 for NH3. 

An additional source of systematic for H 3 + 
can come from the unresolved hyperfinc struc- 
ture (HFS) in combination with possible non- 
th ermal HFS populations in the ISM. As noted 
bv iKeto & Rvbicki I (|2010l) . HFS lines reduce the 
effective optical depth of the molecular rotational 
transition by spreading the emission out over a 
wider bandwidth. To our knowledge, the HFS has 
not been resolved yet for H3 16 + either in labora- 
tory, or astronomical measurements. An expected 
size of the hyperfine splittings can be estimated 
using analogy with ammonia. For the latter the 
main hyperfine splitting is associated with the spin 
of nitrogen I\. The maximum hyperfinc splitting 
caused by the hydrog enic spin I is about 40 kHz 
(|Ho &: Townes 1119831 ). This splitting includes in- 
teraction with the nitrogen spin ~ (Ii -I) and with 
molecular rotation ~ (J ■ I). In hydronium, the 
oxygen nucleus is spinless and there is only inter- 
action with rotation ~ (J ■ I). Hydronium has 
similar electronic structure and close rotational 
constants to ammonia, so its spin-rotational in- 
teraction should be < 40 kHz. Thus, we can 
expect less than 40 m s _1 of the hyperfine band- 
width for 300 GHz lines and smaller splittings at 
higher frequencies. At Tkin ~ 100 K, - a typical 
kinetic temperature of warm and dense gas in the 
star-forming regions, - the thermal width of the 
H 3 + lines is comparable with this HFS splitting. 

For para-H30 + , transitions at 307, 388 and 364 
GHz have 3 HFS components each (two transitions 
with AF — A J are strong, and the remaining one 
is weak). For ortho-H30 + , the 396 GHz transition 
has 9 HFS components (4 strong, 3 weaker, and 2 
weakest). 

The difference in the magnitude of the ener- 
gies of the hyperfine and rotational transitions in 
molecular spectra in the sub-millimiter range is 
considerable: milli-Kelvin for the hyperfine levels 



and tens of Kelvin between rotational levels. This 
means that the HFS levels may be populated ap- 
proximately in statistical equilibrium even if th e 
rotational levels arc not (|Keto fc Rvbicki 1 120101) . 



Therefore, we may suggest that the line centers 
of the sub-millimiter H30 + transitions are not af- 
fected significantly by relative populations of the 
hyperfine levels and that the exepected velocity 
shifts are of a few m s . 

To sum up, the systematic shifts of the line cen- 
ters caused by possible non-thermal HFS popula- 
tions and pressure effects are likely not larger than 
10 m s -1 , which is about 10% of the expected rela- 
tive shift between the para-H30 + Jk = l-i ^ %t 
and 3^ — > 2^ transitions due to ^t-variation. A 
more accurate analysis of the HFS-induced sys- 
tematics will be possible after the HFS is either 
measured, or calculated theoretically. 

Finally, we would like to note that other iso- 
topologues of the hydronium ion also must have 
large sensitivity coefficients Q^. To use them 
we need laboratory studies of the low-frequency 
mixed transitions in the spectra of the partly 
deuterated hydronium ions H2DO + , HD2 + , and 
in D3O+. 

In the near future, high-precision measurements 
in the submillimeter and FIR ranges with greatly 
improved sensitivity will be available with the 
Atacama Large Millimeter/submillimeter Array 
(ALMA), the Stratospheric Observatory For In- 
frared Astronomy (SOFIA), the Cornell Caltcch 
Atacama Telescope (CCAT), and others. Thus, 
any further advances in exploring A^/^t depend 
crucially on accurate laboratory measurements 
(Au/w < 10~ 8 ) of relevant molecular transitions 
in the submillimeter and FIR ranges where reliable 
spectroscopic data are still relatively poor. 

We are grateful to V. Kokoouline for bringing 
inversion spectra of H3 + to our attention. The 
work has been supported in part by the RFBR 
grants 08-02-00460, 09-02-12223, and 09-02-00352, 
by the Federal Agency for Science and Innova- 
tions grant NSh-3769.2010.2, and by the Chinese 
Academy of Sciences, grant No. 2009J2-6. 
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